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Abstract 

We present Gamow- Teller strength distributions from shell model Monte 
Carlo studies of /p-shell nuclei that may play an important role in the pre- 
collapse evolution of supernovae. We then use these strength distributions to 
calculate the electron-capture cross sections and rates in the zero-momentum 
transfer limit. We also discuss the thermal behavior of the cross sections. 
We find large differences in these cross sections and rates when compared 
to the naive single-particle estimates. These differences need to be taken 
into account for improved modeling of the early stages of type II supernova 
evolution. 

PACS numbers: 26.50,+x, 23.40.-s, 21.60Cs, 21.60Ka 
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I. INTRODUCTION 



The impact of nuclear structure on astrophysics has become increasingly important, par- 
ticularly in the fascinating and presently unsolved problem of type-II supernovae explosions. 
The possibility to detect neutrinos ejected before the inf ailing matter reaches the neutrino 
trapping density [1] will allow us, for the first time, to understand whether models of the 
precollapse evolution are in reasonable agreement with observation. These studies will also 
shed light on the effective electron-to-proton ratio in the early stages of the collapse. 

Key inputs for the precollapse scenario are the electron capture cross sections and rates 
for iron group nuclei . The core of a massive star at the end of hydrostatic burning is 
stabilized by electron degeneracy pressure as long as its mass does not exceed the appropriate 
Chandrasekhar mass Mch- If the core mass exceeds Mqh, electrons are captured by nuclei 
@. Thus, the depletion of the electron population due to capture by nuclei is a crucial factor 
determining the initial collapse phase. 

The reduction of the electroweak interaction matrix element to the zero momentum 
transfer limit for the nuclear sector leads directly to the Gamow- Teller operator as a primary 
ingredient in electron-capture cross section calculations. The Gamow- Teller (GT) properties 
of nuclei in the iron region of the periodic table are known to be crucial for supernova physics 
For many of the nuclei that are present in this early stage of the presupernova ||, GT 
transitions contribute significantly to the electron capture cross sections. Due to insufficient 
experimental information, the GT + transitions, which change protons into neutrons, have so 
far been treated only qualitatively in presupernova collapse simulations, assuming the GT + 
strength to reside in a single resonance whose energy relative to the daughter ground state 
has been parametrized phenomenologically Q ; the total GT + strength has been taken from 
the single particle model. Recent (n,p) experiments |5|-|§, however, show that the GT + 
strength is fragmented over many states, while the total strength is significantly quenched 
compared to the single particle model. A recent update of the GT + rates for use in supernova 
simulations assumed a constant quenching factor of 2 || , and included the known low-lying 
transitions in the calculations of the rates. 

In the presupernova collapse, electron capture on p/-shell nuclei proceeds at temper- 
atures between 300 keV and 800 keV 0. While presupernova collapse studies took into 
account the temperature effects on the electron capture rate induced by changes in the 
electron distribution, potential effects on the Gamow- Teller strength distribution have been 
neglected, and use was made of the extreme shifting assumption [H| which states that the 
centroid corresponding to each parent excited state is shifted upward in the daughter nucleus 
by the excitation energy of the parent state. As experiments measure only the ground state 
distribution, this hypothesis has necessarily to be tested theoretically, and it has been found 
approximately valid in a restricted shell model study pT| , [12| ; however, the recent SMMC 
calculations of Ref. clearly indicate changes of the Gamow- Teller strength distribution 
with temperature. Thus it appears desirable to take possible changes of the Gamow- Teller 
strength distribution with temperature also into account when calculating presupernova 
electron capture rates. We will consistently do this in this paper. 

In contrast to previous studies, the calculations presented in this paper solve the full shell 
model problem for the Gamow- Teller strength distributions in the Ohu /p-shell space using 
a realistic residual interaction. We use the shell model Monte Carlo (SMMC) approach to 
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find the Gamow- Teller strength distributions. SMMC has the added advantage that it treats 
nuclear temperature exactly. These distributions are then used to calculate the electron- 
capture cross section as a function of the incident energy of the electron. We also calculate 
the electron-capture rates assuming a relativistic electron gas and a variety of electron 
chemical potentials. In Section II we discuss the methods used in these calculations, while 
in Section III we present our results. A discussion and conclusion follow in Section IV. 



II. THEORETICAL BACKGROUND 

A. Shell Model Monte Carlo techniques 

The SMMC method is based on a statistical formulation of the nuclear many-body 
problem. In the finite-temperature version of this approach, an observable is calculated as 
the canonical expectation value of a corresponding operator A at a given nuclear temperature 
and is given by [0-0 

Tr a [C/] 

where U = exp(—/3H) is the imaginary-time many-body propagator, Tr^U is the canonical 
partition function for A nucleons, H is the shell model Hamiltonian, and (3 = 1/Tjy is the 
inverse temperature. 

The SMMC canonical expectation values are based on the discretization of the many- 
body propagator, exp(—j3H), into a finite number of imaginary time slices, N t , each of 
duration A/3 = (3/N t . At each time slice the many-body propagator is linearized via the 



Hubbard-Stratonovich transformation [p!q , ^9|1 ; observables are thus expressed as path in- 



tegrals of one-body propagators in fluctuating auxiliary fields. The integration is carried 
out by a Metropolis random walk [20]. To circumvent the sign problem encountered in the 
SMMC calculations with realistic interactions, we use the extrapolation procedure outlined 



in Refs. 21.22 



Yet another, but distinct, source of the sign problem is an odd number of nucleons in 
the canonical expectation values [ 17j . We overcome this problem by a number-projection 



technique, first employed in ]23| and subsequently used in JT7|, that allows us to extract 
information concerning odd- A nuclei from the neighboring even-even system. The partition 
function for the odd- A nucleus A' is given by 

Z A , = Jv [a] G{a)Tr A ,U a = Jv[a] G(a)Tr A U a ^^ , (2) 

where we have introduced the trace over the even-even nucleus A. The integration is over 
all auxiliary field variables, a, and G(a) is a Gaussian factor that arises when the Hubbard- 
Stratonovitch transformation is employed. We define our Metropolis sampling weight as 
W(a) = G(a)Tr A U a , and in the case of number projection, the Monte Carlo sign becomes 
$(<r) = Tr A >U a /Tr A U a , which is small, but fairly stable against fluctuations at least down 
to T/v = 0.8 MeV for realistic interactions. We project from the nearest even-even system 
with A' + 1 particles. 
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The SMMC method is also used to calculate the response function -Ra(t) of an operator 
A at an imaginary-time r. The response describes the dynamical behavior of the nucleus 
under the influence of the operator, and contains information about the nuclear spectrum. 
This is seen by using a spectral distribution of initial and final states | i), \ f) with energies 
Ei and Ef 



,2 



= < A = ^~W } = Ei (2J, + l)e~™ 

(3) 

Note that the total strength for the operator is given by R(t = 0). The strength distribution 
is related to Ra( t ) by a Laplace Transform, 

/oo 
S A {E)e~ TE dE. (5) 
-oo 

Note that here E is the energy transfer within the parent nucleus, and that the strength 
distribution S GT +{E) has units of MeV -1 . 



B. Weak interaction cross sections in nuclei 

Electron capture by nuclei is mediated by the electro- weak charged current. At energies 
appropriate for the pre-collapse supernova, the W-boson propagator may be collapsed to 
point coupling and the weak vertex represented by the (V-A) coupling with the universal 
strength g wk . The complexity of computing these channels for electron capture arises from 
the need to treat both the weak sector and the details of the nuclear structure effects with 
suitable finesse to obtain reliable cross sections. 

Starting from the exact expression for the semi-leptonic process corresponding to electron 
capture, one can approximate the three-momentum transfer to the nucleus to be zero in the 
nuclear sector corresponding to computing the nuclear matrix elements between initial and 
final states. This allows us to calculate and predict the nuclear sector with a high level 
of accuracy. The lepton traces and kinematics are still included correctly as is energy 
conservation between the initial and final states. The sacrifice of the magnitude of the 
momentum transfer to the nucleus in the nuclear sector is reasonably well justified at the 
energy domain operative in the precollapse stellar core where energies are low and the energy 
transfer from the incident charged lepton to the nucleus is mainly used for internal nuclear 
excitation. 

Electrons with energy E e interact with the nucleus through the weak interaction. For 
charge exchange interactions 

e-+A(N,Z)^A(N + l,Z-l) + v e -. (6) 
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The parent nucleus absorbs a part of the incident energy and the neutrino carries the re- 
maining energy e v . The energy absorbed by the parent is given by the difference between 
the initial and final nuclear states E = Ef — E iy and the states are labeled by {| /), | i)}, 
respectively. Energy conservation requires that E e = E + e u . The details of the derivation 
of the matrix element are given in p4|. The angular cross section follows as 



da elV 



dQ (2tt) 



E ^Vl £ I if I ^ I j > 

leptonspins 1 Mr Mi 



(7) 



where Hy/ is the weak interaction Hamiltonian, Jj ( Jf) is the initial (final) spin of the nuclear 
state, and Mi, Mf are projections of the angular momentum of the initial and final states 
of the nucleus. In the limit of small three-momentum transfer, the partial cross section for 
electron capture between fixed initial and final states at a given incident electron energy 
reduces to 

= ^ror ' <' " * " <> i 2 < 8) 

where g w k = 1.1661 x 1CT 5 GeV~ 2 is the weak coupling constant, and 

i A 

Am = - 7 =G j ^E t +W°'im(0 • (9) 

Here Ga is the axial vector form-factor which at zero momentum is Ga = 1.25. The Gamow- 
Teller operator is GT + = J2i T t cr ij where is the isospin raising operator, and Oi is the 
Pauli spin operator for nucleon i. The Clebsch-Gordon coefficients that arise when applying 
the Wigner-Eckart theorem give selection rules for this operator that require transitions 
from an initial state with quantum numbers (Jj, Tj, IT) to a final state (J/, Tf, Hf) such that 
AJ = J f - Ji = 0, 1, but no Jf = J i = 0, AT = T f - Ti = 0, 1, and U f = II;. 

Our previous studies showed that the experimentally observed quenching of the total GT 
strength is consistently reproduced by the correlations within the full p/-shell if a renormal- 
ization of the spin operator by the factor 0.8 is invoked p5| , |26[| . The same renormalization 
factor had already been deduced from sd-shell |27fl and p/-shell nuclei with A < 49 |28|j29[ 



and thus appears to be universal. In this work, we also renormalize the spin operator by 
0.8, a eff = 0.8<7i. 

In order to obtain the total cross section at a given incident electron energy, one must 
sum over all initial and final states in Eq. (^|). We are generally interested in situations 
when the initial nucleus is at a finite temperature, so that the initial state sum needs to be 
weighted by the appropriate statistical factors. Thus, the total cross section is 

if Za 

= g& / dm - o 2 e eM ~ m gnaii t^^mm -e,+ . (n) 

if A j = i 

Using the SMMC expression for the Gamow- Teller strength distributions and shifting the 
energy scale in (4) by an appropriate Coulomb correction and the proton-neutron mass 
splitting, the total cross section can be recast into the form 
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°(E e ) = 6gl k J daE e -0 2 j^S GT+ (OF(Z,E e ) . (12) 

Here we have additionally accounted for the distortion of the electron's wave function 
due to the Coulomb field of the nucleus by introducing the Fermi function F(Z, E e ) = 
2nrj/ [exp(2nri) — 1] where r\ = Za//3, a is the fine-structure constant, (3 —\ /3 | is the mag- 
nitude of the velocity of the incoming electron, and Z is the charge of the parent nucleus. 

The presupernova electron capture rate A ec is then given by folding the total cross section 
with the flux of a degenerate relativistic electron gas: 

a. = r^s GT ^T N )-^— rd PP \-t + Ee y F \ z f; ] — - , (is) 

oloosec Jo {m e c z ) Jc 1 + exp [p e {E e — fx e )\ 

where l/f3 e , p, and /i e are the electron temperature, momentum, and chemical potential, 
and L = (Qjj — mlc 4 ) 1 ^ 2 for Qif < —m e c 2 , and otherwise, and Qif = —E = E { — Ef is the 
energy difference between the nuclear states % and /, respectively We have indicated that 
the Gamow- Teller distribution is calculated at the finite nuclear temperature T/v which, in 
principle, is the same as the one of the electron gas T e = 1/ ' (3 e . However, we will study below 
the uncertainties introduced in the electron capture rate due to the temperature dependence 
of the Gamow- Teller distribution by varying T/v for fixed T e . 



III. RESULTS 
A. Validation 



The SMMC calculations of Gamow- Teller strength distributions were detailed in Ref. |T3 
In the present work, we have calculated several more nuclei in addition to those found in 
Ref. |R|. Each even-even nucleus is calculated at a temperature of T N = 0.5 MeV, while 



the ^-extrapolation required to circumvent the SMMC sign problem allows only to cool 
odd-A-nuclei down to a temperature of T N = 0.8 MeV. We use a A(3 = 0.03125 MeV" 1 , 
and approximately four-thousand samples per extrapolation Hamiltonian. We use the KB3 
residual interaction [|30| which is well suited for full Ohu studies throughout the lower pf- 



shell region PH|. Using the numerical techniques described in Ref. [p!| , |rT|l , we construct a 
strength distribution from the response function of the Gamow- Teller operator. Recall from 
Ref. [0 that good agreement was found between experiment and theory for those nuclei 
that had been experimentally investigated. 

In order to demonstrate and validate that our SMMC odd-A projection procedure gives 
reasonable results for the strength distributions at finite temperatures, we compare our 
SMMC results in Sc to standard shell model diagonalization calculations. 45 Sc does not 
actually play a role during the presupernova collapse. Its choice is, however, dictated as 
it is the heaviest odd-A nucleus for which state-by-state calculations of the Gamow- Teller 
strength distribution at finite temperature in the complete p/-shell are feasible. Using Eq. 
3 above, we calculated Sgt+{E,Tn) within the diagonalization approach for all initial states 
in 45 Sc up to an excitation energy of 4.68 MeV. At a nuclear temperature of 0.8 MeV, the 
last excited state considered contributes less than 0.09% to the partition function. Figure 
1 compares the thermal distribution (dotted histogram) obtained by diagonalization with 
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the SMMC result (solid line) calculated at the same temperature T/v = 0.8 MeV. The 
distributions are plotted as function of the nuclear energy transfer E = Ef — E^ required in 
(12). (The Coulomb energy of the final state has been corrected using the shift as defined in 
f25fl .) One observes that the SMMC calculation reproduces both the position and the width 
of the thermal distribution very well. The diagonalization calculation gives a total strength 
of 0.60, while the SMMC result is 0.52±0.08. 

To study the temperature dependence of the Gamow- Teller strength distribution we have 
performed the diagonalization calculations at different temperatures T/v = 0., 0.3 MeV and 
0.8 MeV The results are also given in Fig. 1. At first we note that the T = distribution can 



be compared to the experimental data of Ref. [|52]] . The total calculated strength (0.52±0.08) 
agrees well with the measurement (0.48 up to excitation energies of 9 MeV; the observed 
strength at higher excitation energies is likely due to contributions outside the present model 
space f28|). However, the Gamow- Teller peak of the experimental distribution is about 1 
MeV higher in energy than the calculation, which makes the agreement between data and 



shell model calculation worse than in the other pf -shell nuclei [[26|,|n|. With increasing 



temperature the Gamow- Teller distributions are broadened and its center is shifted slightly 
to lower transfer energies. Both effects have already been noticed in the SMMC studies of 
i"5| . The broadening is mainly due to the fact that the number of states contributing to the 



thermal ensemble increases. The energy shift is caused by the decrease of pairing energy 
in the higher excited states. We mention that the large peak seen in the distributions at 
T/v = and 0.3 MeV corresponds principally to 45 Sc(J = 7/2) — > 45 Ca(J = 5/2) transitions. 
In the distribution at T/v = 0.8 MeV this peak appears lower in energy at E pa 5 MeV and 
its strength is noticeably reduced. 

We also want to understand how the electron capture cross section (11), as calculated 
with the SMMC Gamow- Teller strength distribution at T/v = 0.8 MeV, compares to the 
one obtained from the respective diagonalization results. The middle plot in Fig. 1 shows 
that the SMMC result reproduces the thermal diagonalization data down to E e = 10 MeV 
very closely, and is within 30% of the exact result down to E e = 5 MeV. At very low 
electron energies the cross section becomes very sensitive to the discrete level structure 
of the Gamow- Teller transitions. These weak transitions are not well reproduced by the 
SMMC calculation leading to a noticeable underestimate of the cross section within the 
SMMC approach. However, the cross section is already very low at these electron energies 
and, as we will see next, this shortcoming of the SMMC calculation is not too important 
for the electron capture rate. Note that the standard shell model calculations yield almost 
oscillatory behavior at low E e (dashed line), which reflects its discrete level structure at 
low excitation energies. For completeness, we have also plotted the electron capture cross 
sections calculated for the Gamow- Teller distributions at Tv = and Tv = 0.3 MeV. We 
observe that differences in the distribution only affect the cross section at rather low electron 
energies. 

The quantity of interest for presupernova studies is the electron capture rate (12). In 
the lower panel of Fig. 1 we compare the rates, calculated from the SMMC and from the 
diagonalization Gamow- Teller strength distributions at T/v = 0.8 MeV, as a function of 
electron chemical potential. The agreement is very satisfying for the entire range /i e < 10 
MeV indicating that the differences between the SMMC and diagonalization electron capture 
cross sections at low electron energies are washed out when folded with the electron flux. 
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As mentioned above, SMMC calculations with realistic interactions cannot be performed 
for odd- A nuclei at temperatures below T N = 0.8 MeV. As the relevant temperature during 
the early stage of the presupernova collapse is lower, we have investigated the inaccuracy 
introduced in the electron capture rates if the Gamow- Teller distribution at these lower 
temperatures is approximated by the one at T/v = 0.8 MeV. Choosing the extreme case, 
T/v = 0.3, we have performed one consistent calculation of the electron capture rate using 
the appropriate electron flux and Gamow- Teller strength distribution (determined within the 
diagonalization approach). This result is compared with a calculation in which the Gamow- 
Teller distribution is approximated by the SMMC distribution calculated at T/v = 0.8 MeV. 
We observe a maximal difference of a factor 4 for \x e < 10 MeV which mainly reflects 
deviations in the Gamow- Teller distribution at low energies. Thus, it can be reduced if the 
weak discrete transitions which are experimentally known are added by hand to the SMMC 
results, following the procedure of ||. Obviously the present SMMC approach becomes 
increasingly more accurate for later stages of the collapse where the temperature is higher 
and for even-even parent nuclei where SMMC calculations at the relevant temperature are 
possible. 

Finally, we want to study how much the statistical Monte Carlo uncertainties and the 
errors related to the inversion technique used to determine the Gamow- Teller distribution 
from the response function (5) affect the electron capture rates. The major uncertainty 
introduced in our Gamow- Teller strength distribution is in the position of the centroid which 
has an error of about ±0.5 MeV. In Fig. 2 we show how much such a change would affect 
A ec . Decreasing (increasing) the position of the centroid by 0.5 MeV increases (decreases) 
the rate by an approximate factor of 1.4 (1.3). 



B. Electron capture results 

Ref. fl3| presented SMMC Gamow- Teller distributions for those pf -shell nuclei in the 



iron mass region for which experimental data are available ( 54 ' 56 Fe, 58 > 60 > 62 > 64 Ni ; 51 V, 55 Mn 
and 59 Co). The agreement between theory and data was very satisfying in all cases, and it 
was concluded that the SMMC approach has the predictive power to estimate the Gamow- 
Teller strength distribution for nuclei for which no data exist. We have extended the study 
of Ref. |I3| by performing SMMC calculations for several of those p/-shell nuclei which are 



of interest for the presupernova collapse: 55,57 Co, 56 Ni, 50 ' 52 Cr, 55 > 58 Fe, 56 Ni, and 50 Ti. The 
Gamow- Teller distributions for these nuclei are plotted in Fig. 3 as a function of energy 
transfer E = Ef — Ei. As for 45 Sc, the Coulomb correction has been performed following the 
prescription given in Ref. pp| . Each Gamow- Teller distribution has been interpolated onto 



a fine energy grid with a spline function. The spline function was then used to calculate the 
cross section using the formulae described in Section II. 

Before presenting results for the electron capture cross sections and rates, we note that 
our approach, as presented in Section II, implicitly assumes that the calculation reproduces 
the mass splitting between parent and daughter nucleus correctly. As has been shown in 
|26| , |25[| the KB3 interaction indeed describes masses of p/-shell nuclei reasonably well (typi- 
cally within 500 keV of the data; as the deviations have a systematical trend the inaccuracy 
on the relevant mass splitting between parent and daughter should be even smaller). Many 
of the SMMC masses presented in [EB1 have very recently also been calculated within the 
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conventional diagonalization approach, and the SMMC results (after subtraction of the ap- 
propriate internal excitation energy) have been confirmed in all cases [P3| . These studies, 



however, clearly showed that the KB3 interaction overbinds the /V/2 subshell closure. 

This shortcoming can be circumvented [34j by using the FDP6 interaction |35 ]. However, 
this force does not reproduce the Gamow- Teller strength correctly ]36|]. Weighing the advan- 



tages and disadvantages of the two interactions, we have chosen the KB3 interaction with 
the price that we had to shift the energy scale of the Gamow- Teller distributions in order 
to correct for the f'7/2 subshell overbinding. As can be seen in Ref. Jl3[| , after correction 
the SMMC calculation reproduces the observed Gamow- Teller distributions in 54 Fe and the 
various Ni isotopes very well. 

In Table 1 we compare the SMMC results for the Gamow- Teller peak in the daughter 
nucleus with the parametrization given by Fuller, Fowler, and Newman Q and subsequently 
used in the update of the presupernova electron capture rates by Aufderheide et al. ||. 
We observe that Ref. [|[] generally places the Gamow- Teller strength for even-even parent 
nuclei at too high an excitation energy in the daughter, while it is the opposite for odd- /I 
parent nuclei. The same trend as in the SMMC distributions is found in the data and has 
already been pointed out in Ref. p7 |. From the SMMC calculations and the data, we can 



conclude that for an even-even parent, the Gamow- Teller strength, at low temperatures, is 
at low daughter excitation energies (~ 2 MeV), while it is at significantly higher excitation 
energies (~ 5 MeV) for an odd-A nucleus. This different behavior is related to the pairing 
energy difference between the even-even parent and odd-odd daughter. As the Gamow- Teller 
distributions usually have widths of order 1-2 MeV, our SMMC calculations for odd-A nuclei 
might miss weak Gamow- Teller transitions at low excitation energies which could potentially 
increase the electron capture rates. It appears therefore reasonable to follow the prescription 
of [Q and add experimentally known transitions at low energies to the rates obtained from 
the SMMC Gamow- Teller distributions for odd-A nuclei. Such a procedure seems to be 
unnecessary for even-even parent nuclei. 

The electron capture cross sections calculated from the Gamow- Teller matrix elements 
are presented in Fig. 3 as a function of the incident electron energy E e . The general behavior 
of the cross section can be easily understood. To trigger the electron capture process requires 
a minimum electron energy given by the mass splitting between parent and daughter, Qif. 
(This threshold is lowered by the internal excitation energy at finite temperature.) In even- 
even parent nuclei the Gamow- Teller strength, at low temperatures, is centered at daughter 
excitation energies of order 2 MeV. Thus the electron capture cross sections for these parent 
nuclei increase drastically within the first couple of MeV of electron energies above threshold, 
reflecting the Gamow- Teller distribution. With increasing electron energies it continues to 
raise modestly caused by the (E e — £) 2 factor in the cross section integral (12). As Qif 
increases with neutron excess, caused by the increase of the symmetry energy, electron 
capture cross sections, at fixed electron energies, decrease with neutron excess. In odd- 
A nuclei the Gamow- Teller distribution, at low temperatures, peaks at noticeably higher 
daughter excitation energies. Thus the electron capture cross sections are shifted to higher 
electron energies for odd-A nuclei in comparison to even-even parent nuclei by about 3 MeV, 
reflecting the difference in the Gamow- Teller peak positions. 

Shown in Fig. 5 are the electron capture rates as a function of the electron chemical 
potential. We have used T e = T N = 0.5 MeV for the electron capture rates on even parent 
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nuclei, while for odd- A parent nuclei we used T/v = 0.8 MeV for the reasons discussed above. 
The electron chemical potential depends on the electron density, the electron fraction in the 
medium, and the temperature of the electrons. In the precollapse phase of supernova, a 
reasonable approximation is given by f38 



/i e = l.ll(p 7 F e ) 1/3 



7T \ 2 Tl 



111^ (P 7 ^e) 2/3 



-1/3 

(14) 



where T e is in units of MeV, and the electron density p-j = p/10 7 . Thus, for the nuclei studied 
here p e ~ 1.5 — 2 MeV under typical presupernova conditions (T/v ~ 0.4 MeV, pj ~ 6). 

Do the present electron capture rates indicate potential implications for the pre-collapse 
evolution of a type II supernova? To make a judgement on this important question, we 
compare in Table II the SMMC rates for selected nuclei with those currently used in collapse 
calculations ||. For the comparison we choose the same physical conditions as assumed in 
Tables 4-6 in [BJ. Furthermore the table lists the partial electron capture rate which has 
been attributed to Gamow- Teller transitions in Ref . || . At first we note that for even parent 
nuclei, the present rate approximately agrees with the currently recommended total rate. A 
closer inspection, however, shows significant differences between the present rate and the 
one attributed to the Gamow- Teller transition in ||. As discussed above, the origin of this 
discrepancy is due to the fact that Ref. [|J places the Gamow- Teller resonance for even-even 
nuclei systematically at too high an excitation energy. Of course, this shortcoming has been 
corrected for in Refs. [|],[3| by adding an experimentally known low-lying strength on top 
of the one attributed to Gamow- Teller transitions. However, the overall good agreement 
between the SMMC results for even-even nuclei and the recommended rates indicates that 
our SMMC approach also accounts correctly for this low-lying strength. This has already 
been deduced from the good agreement between SMMC Gamow- Teller distributions and 
data including the low-energy regime. We conclude therefore that, for even-even nuclei, the 
SMMC approach is able to predict the total electron capture rate rather reliably, even if 
no experimental data are available. We note that our SMMC rate is somewhat larger than 
the recommended rate for 56 Fe and 60 Ni. In both cases the experimental Gamow- Teller 
distribution is known and agrees well with the SMMC results [|13j. While the proposed 
increase of the rate for 60 Ni is not expected to have noticeable influence on the pre-collapse 
evolution, the increased rate for 56 Fe makes this nucleus an important contributor in the 
change of Y e during the collapse (see Table 15 of 0). 

Turning now to electron capture on odd- A nuclei. Here we observe that the SMMC elec- 
tron capture rate, derived from the Gamow- Teller distributions, is significantly smaller than 
the recommended total rate. As already stressed above, this is simply due to the fact that 
for odd- A nuclei the Gamow- Teller transition peaks at rather high excitation energies in the 
daughter nucleus. The electron capture rate on odd- A nuclei is therefore carried by weak 
transitions at low excitation energies. Comparing our rates to those attributed to Gamow- 
Teller transitions in Refs. reveals that the latter have been, in general, significantly 
overestimated which is simply caused by the fact that the position of the Gamow- Teller 
resonance is usually put at too low excitation energies in the daughter (see Table I). For ex- 
ample, Ref. M attributes 87% of the total capture rate on 55 Co to Gamow- Teller transitions 
(upper part of Table II), while our calculation predicts the contribution of the Gamow- Teller 
strength distribution to be negligible. Thus, we recommend that the capture rate on 55 Co 
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is significantly smaller (only 13% of the rate given in Table 4 of Ref. ||). Our SMMC cal- 
culation implies that the Gamow- Teller transitions should not contribute noticeably to the 
electron capture rates on odd- A nuclei at the low temperatures studied in Tables 14-16 in 
@. Thus, the rates for odd- A nuclei given in these tables should generally be replaced by 
the non-Gamow- Teller fraction. 



IV. DISCUSSION 

In this paper we presented a detailed description of the electron-capture cross section and 
rates calculated from the Gamow- Teller distributions obtained using the SMMC method. 
We demonstrated the validity of our odd- A sampling technique in 45 Sc, and we found very 
good agreement between the SMMC and the traditional approach for solutions of the shell 
model. Furthermore, the SMMC approach reproduces the measured Gamow- Teller strength 
distributions very well. In accord with the data, we find an odd-even dependence of the 
Gamow- Teller peak position in the daughter nucleus: while it is generally at low daughter 
excitation energy for even-even parent nuclei (~ 2 MeV), it is at noticeable higher energies 
for odd- A parents (~ 5 MeV). These systematics are not reproduced by the parametrization 
of the Gamow- Teller resonance as adopted for the derivation of the currently recommended 
rates 0,0). This parametrization places the Gamow- Teller resonance usually at too high 
excitation energies for even-even parent while the position is too low for odd- A parents. 
However, this shortcoming has been mainly overcome in the recommended rates by adding 
experimentally known strength at low excitation energies to the Gamow- Teller strength. 

The presupernova collapse occurs at finite temperature and our SMMC approach, for 
the first time, allows to take thermal effects consistently into account. With increasing 
temperature, the Gamow- Teller distribution is broadened and shifted to lower transfer en- 
ergies. However, at the rather low temperatures at which electron capture occurs on nuclei 
in the iron mass region, the temperature dependence of the Gamow- Teller strength distri- 
bution is rather unimportant. For even-even nuclei the distribution does not change too 
much at T N < 0.6 MeV due to the large pairing gap. For odd- A nuclei, the Gamow- Teller 
strength resides at such high excitation energies that the modest temperature dependence 
of the strength does not affect the total rate. Thus it also does not matter that numerical 
problems do not allow us to cool odd- A nuclei below T/v = 0.8 MeV. However, we like to 
stress that the temperature dependence of the Gamow- Teller strength will become impor- 
tant at later stages of the presupernova collapse when temperatures of order T N = 1 MeV or 
higher are reached, as then the strength for both even-even and odd-A nuclei is noticeably 
shifted downwards in transfer energy. At these temperatures the Y e value, however, has 
decreased enough so that electrons are captured on nuclei with Z < 40 and N > 40 for 
which the Gamow- Teller transitions are naively expected to be Pauli-blocked. However, this 
Pauli-blo eking can be overcome by thermal excitation |$9| and by pairing [JD], which, at low 
temperatures, is the more important effect. 

As the Gamow- Teller strength resides at rather high excitation energies for odd-A par- 
ents, the Gamow- Teller contribution is generally negligible in the total electron capture 
rate. This finding is at variance with the recent compilations which propose noticeable 
Gamow- Teller fractions for the capture rate on odd-A nuclei like 55,59 Co. We recommend to 
use only the non-Gamow- Teller fraction of the compiled rates for odd-A nuclei. For even- 
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even parents, the Gamow- Teller distribution is located at very low excitation energies and, 
within our approach, should account for the total electron capture rate. In fact, our rate 
agrees reasonably well with the compiled total rates. However, for nuclei like 56 Fe and 60 Ni, 



the Gamow- Teller strength resides, in agreement between data and SMMC calculation [13 



at such low energies that the present rate is significantly larger than the one given in the 
compilations. 

In summary, the present SMMC approach allows for the first time a microscopic de- 
termination of the Gamow- Teller contributions to the presupernova electron capture rates. 
The present calculation will be extended to other nuclei of potential importance during the 
collapse phase of type II supernova. The SMMC approach can also be extended to heavier, 
more neutron-rich nuclei which will carry the electron capture process at later stages. For 
these nuclei we expect the ability of our method to consistently account for finite tempera- 
ture effects to be quite important. 
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TABLES 



TABLE I. Comparisons of the positions of the Gamow- Teller peaks as calculated within the 
SMMC approach with the parametrization of Ref. Q. The proton and neutron numbers refer to 
the parent nucleus, while Eqt is the position of the Gamow- Teller peak in the daughter nucleus. 
(To obtain this number we have shifted the transferred energy E by the experimental mass splitting 
between parent and daughter nucleus.) £ffn denotes the position of the Gamow- Teller resonance 
due to the parametrization given in Ref. Q and subsequently used in presupernova collapse studies. 
The SMMC calculations have been performed at T/v = 0.5 MeV for even nuclei and at Tjy = 0.8 
MeV for odd- A nuclei. For 50 Cr the single-particle model allows transitions into the fan and f§M 
orbitals. 



N 


Z 


Egt 




28 


28 


2.6 


3.78 


30 


28 


2.0 


3.76 


32 


28 


1.0 


2.0 


28 


27 


6.9 


3.78 


30 


27 


4.8 


3.77 


32 


27 


4.2 


2.0 


28 


26 


2.8 


3.80 


29 


26 


6.1 


5.4 


30 


26 


1.5 


3.78 


32 


26 


-0.5 


2.0 


30 


25 


4.7 


3.79 


26 


24 


2.5 


2.0,8.7 


28 


24 


1.3 


3.82 


28 


23 


5.5 


3.83 
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TABLE II. Comparisons of the present SMMC electron capture rates with the total (A ec ) and 
partial Gamow- Teller (A^. T ) rates as given in Ref. ||. Physical conditions at which the comparisons 
were made are given in the last column. 



nucleus 


A cc (sec -1 ) 


A ec (sec 1 ) 


A e G c T (sec" 1 ) 




conditions 








(SMMC) 


(Ref. |) 


(Ref. g) 












55 Co 


2.25E-04 


1.41E-01 


1.23E-01 


Pi 


= 5.86, T 9 


= 3.40 


Ye 


= 0.47 


57 Co 


1.97E-06 


3.50E-03 


1.31E-04 


Pi 


= 5.86, T 9 


= 3.40 


Ye 


= 0.47 


54 Fe 


4.64E-05 


3.11E-04 


9.54E-07 


Pi 


= 5.86, T 9 


= 3.40 


Y e 


= 0.47 


55 Fe 


7.22E-09 


1.61E-03 


1.16E-07 


Pi 


= 5.86, T 9 


= 3.40 


Y e 


= 0.47 


56 Ni 


1.96E-02 


1.60E-02 


6.34E-03 


Pi 


= 5.86, T 9 


= 3.40 


Y e 


= 0.47 


58 Ni 


6.02E-04 


6.36E-04 


4.04E-06 


Pi 


= 5.86, T 9 


= 3.40 


Y e 


= 0.47 


60 Ni 


1.34E-05 


1.49E-06 


4.86E-07 


Pi 


= 5.86, T 9 


= 3.40 


Y e 


= 0.47 


59 p„ 


9 D^F D7 


9 DQF 04 


U.O 1 Ej-UO 


Pi 


= 10.7, r 9 


= 3.65, 


y 
J-e 


— n 


57 Co 


1.21E-05 


7.65E-03 


3.69E-04 


Pi z 


= 10..7, T 9 


= 3.65, 


Y e 


= 0.455 


55 Fe 


6.48E-08 


3.80E-03 


5.51E-07 


Pi 


= 10.7, T 9 


= 3.65, 


Y e 


= 0.455 


56 Fe 


5.86E-06 


4.68E-07 


6.60E-10 


Pi 


= 10.7, T 9 


= 3.65, 


Y e 


= 0.455 


54 Fe 


2.26E-04 


9.50E-04 


3.85E-06 


Pi 


= 10.7, T 9 


= 3.65, 


Y e 


= 0.455 


51y 


6.17E-07 


1.24E-05 


9.46E-09 


Pi 


= 10.7, T 9 


= 3.65, 


Y e 


= 0.455 


52 Cr 


5.22E-07 


2.01E-07 


1.59E-10 


Pi 


= 10.7, T 9 


= 3.65, 


Y e 


= 0.455 


60 M 


7.48E-05 


7.64E-06 


2.12E-06 


Pi 


= 10.7, T 9 


= 3.65, 


Y e 


= 0.455 
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FIGURES 



FIG. 1. Top panel: The Gamow- Teller transitions as a function of the energy transfer 
E = Ef — Ei using SMMC for 45 Sc, and comparing to diagonalization calculations (DD) at 
Tjy = 0.33 and 0.8 MeV. Middle panel: The cross section calculated from both SMMC and diago- 
nalization. Bottom panel: Calculated rates from SMMC (T\r = 0.8 MeV), assuming T e = 0.3 and 
0.8 MeV, compared to the exact results (T e = Tjv). 

FIG. 2. The error band in the rates calculations for 45 Sc, where the position of the calculated 
Gamow- Teller centroid has been shifted by 0.5 MeV. 

FIG. 3. Shown are the GT + distributions for various nuclei studied here as a function of the 
energy transfer from the electron, E = Ef — Ei. These curves were obtained after MaxEnt recon- 
struction of the Gamow- Teller response functions that were calculated in the SMMC framework. 

FIG. 4. Shown are the electron capture cross sections as a function of the incident electron 
energy E e for all nuclei in this study. 

FIG. 5. Shown are the electron capture rates as a function of the electron chemical potential 
energy \i e for all nuclei in this study. All electron temperatures were fixed at T e = 0.5 MeV. 
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